scnPtr = img + (SA)f sM+S . yWinM- I ) *pixPerScn ^fe . xOf sM ; 
if (levels==8) { W W 

if (S . segmntMode==3) { ./* directed acuity mode 3; pack pixels */ 
/+ S.xWinM and S.yWinM known to be mod 8 Sc. mod 2, respectively */ 
for (j=0; j<S.yWinM; j+=2, scnPtr-=2*pixPerScn) { 
pxlPtr = scnPtr; 

for (i=0; i<S.xWihM; i+=2, pxlPtr+=2) { 
bA = (int) *pxlPtr; 
bB = (int) * (pxlPtr+1) ; 
bC = (int) * (pxlPtr-pixPerScn) ; 
bD = (int) * (pxlPtr-pixPerScn+1) ; 

/* bit 0 is set if halftone (&1) */ 
/* bit 1 is set if verticle (&2) */ 

/* code for diagionals is same value in 224 and 28 Iocs */ 
/* diagionals made only when bit d is blacker than a,b,c */ 
vector = 0 ; 

vector =8; /* halftone A */ 
vector =4; /* halftone B */ 
vector =2; /* halftone C */ 
vector =1; /* halftone D */ 
switch (vector) { 

case 0:if/* ( ( (bD&224) < (bB&224) ) && ( (bD&224 ) < (bC&224 ) ) ) */(0 
byt = (bD&227) | ( (bD>>3) &28) ; break; 



if 
if 
if 
if 



(bA&l) 
(bB&l) 
(bC&l) 
(bD&l) 



) 1 



else 
byt 





} 


byt 




case 


1 : 


byt 








byt 




case 


2 : 


byt 








byt 




case 


3 : 


byt 








byt 




case 


4 : 


byt 








byt 




case 


5 : 


byt 








byt 




case 


6: 


byt 




case 


7: 


byt 








byt 




case 


8:if/*( ( (1 




} 


byt 






else 


i" 




byt 






} 


byt 




case 


9: 


byt 








byt 




case 


10: 






case 


11 : 


byt 








byt 




case 


12 : 






case 


13 : 


byt 





(bA&2)? (bA&227) | ( (bB>>3) &28) 

: (bA&227) | ( (bC>>3) &28) ; 
grad(byt); break; 

(bA&2)? (bA&227) I ( (bB>>3)&28) 

: (bA&227) I ( (bC>>3) &28) ; 
grad(byt); break; 

(bB&2)? (bA&224) I ((bB>>3)&28) | (bB&3) 

: (bA&224) |m28 (bA) | (bB&3) ; 
grad(byt) ; break; 
(bA&2)? (bA&227) I ( (bB>>3)&28) 

: (bA&227) |m28 (bA) ; 
grad(byt) ; break; 
(bC&2)? (bA&224) Im28 (bA) | (bC&3) 

: (bA&224) I ( (bC>>3) &28) I (bC&3) ; 
grad(byt) ; .break; 
(bA&2)? (bA&227) Im28 (bA) 

: (bA&227) I ( (bC>>3) &28) ; 
grad(byt) ; break; 

bA&227 Im28 (bA) ; break; 
bA&227 |m28 (bA) ; 
grad(byt); break; 



(bD&227) | ( (bD>>3) &28) ; break; 

(bD&2) ? m224 (bC) I ( (bB>>3) &28) I (bD&3) 

: m224 (bB) | ( (bC>>3) &28) | (bDS:3) ; 
grad(byt); break; 

(bB&2) ? m224 (bC) I ( (bB>>3) &28) I (bB&3) 

: m224 (bB) I ( (bC>>3) &28) | (bB&3) ; 
grad(byt); break; 

m224 (bB) | ( (bB>>3)&28) | (2) ; 
grad (byt) ; break; 

m224 (bC) | ( (bC>>3) &28) ; 
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^fcbyt = grad(byt); break; 
case 14^jyt = (bD&227) | ( U >3)&28); break; 

case 15: byt = bA; break; 

tempi = (uch)byt; 

putc (tempi, stmPtr) ; /* even is lineart */ 

} else { /* mode must be 0,1 or 2 */ 
/* S.xWinM known to be mod 4 */ 
for (j=0; j<S.yWinM; scnPtr- =pixPerScn) { 

pxlPtr = scnPtr; 

for (i-0; i<S.xWinM; i++, pxlPtr++) { 
tempi =' *pxlPtr; 

putc(templ, stmPtr); /* even is lineart */ 

} else { /* dont do segmentation for quants less than 8 */ 
/* quantize the value, then map into a byte, 0 - 255 */ 
qdblO = (dbl) (l«levels) ; 
qdbll = (qdblO - .000001) / 255.0; 
qdbl2 = 255.0 / (qdbl0-1.0); 
/* S.xWinM known to be mod 4 */ 
for (j-0; j<S.yWinM; scnPtr-=pixPerScn) { 

pxlPtr = scnPtr; 

for (i=0; i<S.xWinM; i++, pxlPtr++) { 
tempi = *pxlPtr; 
if ( (templ==255) | | (tempi— 0) ) { 

putc(templ, stmPtr); 
} else { 

qdbl3 = ( dbl) tempi; 

qdbl4 = ( floor ( qdbl3 *qdbll )) *qdbl2; 
\ P utc( ( (unsigned char) qdbl4) , stmPtr); 

\ close (stmPtr) ; 

} 



it 

rad(byte) 
at byte; 
return byte; } 



radl (byte) 

a /* b re turns byte with entrys not equal; ie a gradient */ 
int bl, b2, b3, b4, b5, vector; 

bl = byte&224; b2 = (byte&28 ) <<3 ; 
if ( ! (bl==b2) ) { 

return byte; 
} else { 

vector = bl>>5; 
switch (vector) { 

case 7: b3=24; break; /* 224 */ 
case 6: b3=28; break; /* 192 */ 
case 5: b3=24; break; /* 160 */ 
case 4: b3=20; break; /* 128 */ 
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case 3: b3=8 

case 2: b3=4 

case 1: b3=0 

case 0: b3=4 



^ak; /* 96 */ 

Wiak; /* 64 */ 

break; /* 32 */ 

break; /* 0 */ 



b4 = byte&3; 
b5 = bl|b3|b4; 
return b5; 



28 (byte) 
nt byte; 

if (byteSd28) return 28; else return 0 



-24 (byte) 
nt byte; 

if (byte&128) return 224; else return 



.sassemr » 



srcPtr = srcStrt; 
ofsPtr = srcStrt + xwin*2 
nxtPtr 

IstPtr ■ srcStrt - xwin*4 
for (i=0; i<ywin; i++) { 
for (j=0; j<xwin; j++) 



/* offset by one raster */ 
srcStrt + xwin*4; /* offset by two rasters */ 
/* offset by two rasters */ 



byteP = *srcPtr; 
byteQ - *(srcPtr+2); 
byteR = *nxtPtr; 
byteS - MnxtPtr+2); 

byteU = *lstPtr; /* U & V have been unpacked last raster */ 
bytev = * (srcPtr-2) 
crnrP = * (lstPtr-2) 
crnrQ = * (lstPtr+2) 
crnrR = * (nxtPtr-2) 
crnrS = *(nxtPtr+2) 
dglPS = ( (crnrP&l) &&(crnrS&l) ) ; /* 
dglQR = ( (crnrQ&l) && (crnrR&l) ) ; /* 



both PS- corners lineart 
both QR corners lineart 



/* lsb indicates lineart; next lsb indicates horizontal 
if (byteP&l) { /* p is lineart - gen pqrs as lineart pixels */ 
tmpp = (byteP&224) ; 
pxlp = tmpp | (tmpp>>3) ; 

if ( ( (byteQ&2) " (byteR&2) ) /* one vert, the other horz */ 

&&{byteQ&l) &&(byteR&l) /* and both sides are lineart * 

&&(byteU&l) &&(byteV&l) /* and four sides are lineart * 

&&{(dglPS) || (dglQR) ) ) { /* and either ps or qr lineart 

/* if both diags lineart, replicate lowest along diag. */ 

/* if only one diag lineart, replicate in dir of diag. */ 

pxla = pxlp | 1 ; 

tbot = byteP&28; 

pbot = tbot | (tbot<<3) | 1; 

if ( (dglPS) &&(dglQR) ) { /* both diagionals lineart */ 
/* replicate the lowest */ 

if (pxla<pbot) { /* replicate a to d; b or c interpolated 
pxld = pxla; 

if (byteP&2) { /* p horz; bot of p goes to c; b intrp */ 
tmpq ■ (byteQ&224) ; ' 
pxlq = tmpq | (tmpq>>3) ; 
pxlb - { (pxlp+pxlq)/2) |1; 
pxlc = pbot; 

} else { /* p vert; bot of p goes to b; c intrp */ 
tmpr » <byteR&224); 
pxlr * tmpr | (tmpr>>3) ; 
pxlc = ( (pxlp+pxlr) /2) |i; 
pbot; 



pxlb 



} else { /* replicate b to c; d interpolated */ 

pxlb = pbot; 

pxlc » pbot; 

tmpq = (byteQ&224) ; 

pxlq = tmpq| (tmpq>>3 > ; 

tmpr = (byteR&224); 

pxlr = tmpr | (tmpr>>3) ; 

• pxld = ( (pxlq+pxlr)/2) |l; 

} else { /* only one diagional PS or QR is lineart */ 
if (dglPS) { /* only diagional PS is lineart */ 
/* replicate along PS direction */ 
pxld « pxla; 

if (byteP&2) { /* p is horz */ 
pxlc = pbot; 



©1998 Xerox Corporation 



pJR = cmpq i icmpq>>3) ; 
pxlb = ( (pxlp+pxlq) /2) | 1; 
} else { /* p is vert */ 
pxlb = pbot ; 
tmpr = (byteR&224) ;. 
pxlr = tmpr | (trnpr>>3) ; 
pxlc = ( (pxlp+pxlr) /2) | 1; 

} else { /* only diagional QR- is lineart */ 
/* replicate along QR direction */ 

pxlb = pbot / 

pxlc = pbot / 

tmpq = (byteQ&224) ; 

pxlq = tmpq | ( tmpq>>3) ; 

tmpr = (byteR&224) ; 

pxlr = tmpr | (tmpr>>3); 

pxld = ( (pxlq+pxlr) /2) |1; 



i = i; /* breakpoint */ 
} else { 

pxla = pxlp|l; 

if (byteP&2) { /* p is horz */ 
tmpc = byteP&28; 
pxlc = tmpc | (tmpc<<3) |l; 
if (byteQ&l) { /* q is lineart */ 
tmpq = (byteQ&224); 
pxlq = Cmpq| (tmpq>>3) ; 
if (byteQ&2) { /* q is horz */ 
tmqs = byteQ&28; 
pxqs = tmqs | (tmqs<<3 ) ; 
} else { /* q is verticle */ 

if (byteS&l) { /* s is lineart */ 
tmps » (byteS4224) ; 
pxls = tmps j (tmps>>3) ; 
pxqs = (pxlq+pxls) /2 ; 
} else { /* s is halftone */ 
pxqs = pxlq; 

} } 

pxlb = ( (pxlp+pxlq) /2) I 1; 
pxld » ( (pxlc+pxqs) /2) 1 1; 
} else { /* q is halftone */ 
pxlb » pxla; 
pxld « pxlc; 

} else { /* p is vert */ 
tmpb * byteP&28; 
pxlb - tmpb | (tmpb<<3) |l; 
if (byteR&l) { /* r is lineart */ 
tmpr - (byteR&224) ; 
pxlr * tmpr | (tmpr>>3) ; 
if (byteR&2) { /* r is horz */ 

if (byteS&l) { /* s is lineart */ 
tmps = (byteS&224) ; 
pxls = tmps | (tmps>>3) ; 
pxrs = (pxlr+pxls) /2 ; 
} else { /* s is halftone */ 
pxrs = pxlr; 

} else { /* r is verticle */ 
tmrs = byteR&2 8; 
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pxlc = ( (pxlp+pxlr) /2) j 1; 
pxld = ( (pxlb+pxrs) /2) 1 1; 
else { /* r is halftone */ 
pxlc = pxla; 
pxld = pxlb; 



/* p is halftone 

(byteQ&l) ? pxlp : 
(byteR&l) ? pxlp : 
(byteS&l) ? pxlp : 

pxlp ; 
( (pxlp+pxlq) /2) ; 
( (pxlp+pxlr) /2) ; 



■ generate pqrs as halftone pixels */ 
byteP&254; 

byteQ&254; /* q lineart? */ 
byteR&254; /* r lineart? */ 
byteS&254; /* s lineart? */ 



} else 
pxlp 
pxlq 
pxlr 
pxls 
pxla 
pxlb 
pxlc 

pxld ■ ( (pxlp+pxlq+pxlr+pxls) /4) ; 

*(srcPtr+0) = (uch)pxla; 
*(srcPtr+l) = (uch)pxlb; 
*(ofsPtr+0) = (uch)pxlc; 
*(ofsPtr+l) = (uch)pxld; 

srcP.tr += 2; /* advance by one source pixel */ 
ofsPtr +~ 2; /* advance by one source pixel */ 
nxtPtr += 2; /* advance by one source pixel */ 
lstPtr +« 2; /* advance by one source pixel */ 

} 

srcPtr += xwin*2; /* advance by one more source raster */ 
ofsPtr += xwin*2; /* advance by one more source raster */ 
nxtPtr += xwin*2; /* advance by one more source raster */ 
lstPtr += xwin*2; /* advance by one more source raster */ 
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